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ABSTRACT 

Aims. For more than 1.5 years we spectroscopically monitored the star KELT-6 (BD+31 2447), which is known to host the transiting 
hot-Saturn KELT-6 b, because a previously observed long-term trend in radial velocity time series suggested that there is an outer 
companion. 

Methods. We collected a total of 93 new spectra with the HARPS-N and TRES spectrographs. A spectroscopic transit of KELT-6 b 
was observed with HARPS-N, and simultaneous photometry was obtained with the IAC-80 telescope. 

Results. We proved the existence of an outer planet with a mininum mass M p sin /=3.71±0.21 Mj up and a moderately eccentric 
orbit (e = 0.21ffjj5?) of period P~3.5 years. We improved the orbital solution of KELT-6b and obtained the first measurement of 
the Rossiter-McLaughlin effect, showing that the planet has a likely circular, prograde, and slightly misaligned orbit with a projected 
spin-orbit angle of J=-36± 11 degrees. We improved the KELT-6 b transit ephemeris from photometry and provide new measurements 
of the stellar parameters. KELT-6 appears as an interesting case for studying the formation and evolution of multi-planet systems. 
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1. Introduction 


Analysing stellar radial velocity (RV) time series is an effective 
method of detecting and characterizing distant planets with or¬ 
bital periods of a few years, despite the large observing time 
span required. To look for such companions, particularly inter¬ 
esting targets are stars with transiting planets, because they are 
an ideal laboratory for studying the architecture of multi-planet 
systems. In fact, the orbital geometry of the transiting planet can 
be described through the Rossiter-McLaughlin (RM) effect (e.g. 
Ohta et al. l2005i) by measuring the projected spin-orbit angle and 
spotting the direction of the motion with respect to that of the 
stellar rotation. Up to now, long-term trends have been observed 
in the RVs of a larg e samp le of stars with and without evidence 
of turnover (e.g. lKnutson et ah! 201 4). Among them is KELT-6, 
a late F-type, metal-poor star (V=10.3 mag) hosting the tran¬ 
siting Saturn-mass planet KELT-6 b discovered by the KELT- 


* Based on observations made with (;) the HARPS-N spectrograph on 
the Italian Telescopio Nazionale Galileo (TNG), operated on the island 
of La Palma by the INAF - Fundacion Galileo Galilei (Spanish Observa¬ 
tory of Roque de los Muchachos of the IAC); (it) the Tillinghast Reflec¬ 
tor Echelle Spectrograph (TRES) on the 1.5-meter Tillinghast telescope, 
located at the Smithsonian Astrophysical Observatory’s Fred L. Whip¬ 
ple Observatory on Mt. Hopkins in Arizona; (Hi) the IAC-80 telescope 
at the Teide Observatory (Instituto de Astroffsica de Canarias, IAC). 
Tables 2 and 3 are made available on-line at the CDS (http://cdsweb.u- 
strasbg.fr/). 


North survey (ICollins et al.ll20l4 hereafter Col4). Col4 were 
also able to observe an unexplained residual trend in the RVs 
over a limited time span of 475 days. To understand the cause 
of this acceleration, a spectroscopic follow-up was carried out in 
the framework of the Global Architectures of Planetary Systems 
(GAPS) projecfl using the HARPS- N spec trograph (resolving 
power R = 115 000: ICosentino et al.ll2012l) . We also collected 
new RV data with the Tillinghast Reflector Echelle Spectrograph 
(TRES) spectrograph (R = 44 000; lFureszl l 2008 ). extending the 
total observing time span to 1178 days. 

Together with new photometric measurements, we present 
results that noticeably extend knowledge about the KELT-6 sys¬ 
tem. 


2. Observations and data reduction methods 

We collected 71 HARPS-N spectra (exposure 900 s, typical 
signal-to-noise per pixel S/N~60 at 5500 A) between 2014 
February 9 and 2015 July 3, 31 of which were obtained on 
2015 April 11 during a transit of KELT-6 b and used to study 
the RM effect. The Th-Ar simultaneous calibration was not used 
to avoid contamination by the lamp lines. The spectra and the 
RV measurements were reduced using the latest version (Nov. 
2013) of the HARPS-N instrument Data Reduction Software 

1 http://www.oact.inaf.it/exoit/EXO-IT/Projects/Entries/201 l/12/27_GAPS.html 
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pipeline and applying a G2 mask. The measurement of the RVs is 
based on the weighted cross-correlatio n function (CCF) method 
dBaranne et alJl!996t iPene et aDl 2002 l) . 

With TRES we collected 22 spectra between 2013 December 
13 and 2015 May 27. They were extracted following the pro¬ 
cedures described by iBuchhave et al] d2010l) . The relative RVs 
were derived by cross-correlating the spectra against the highest 
S/R spectrum in the wavelength range 4050- 5650 A. 

Simultaneously with the RM effect measurements gathered 
with HARPS-N, we collected the transit light curve with the 
IAC-80 0.82-m telescope. Data were taken from 21:18 UT to 
5:32 UT, using the CAMELOT camera (E2V 2kx2k CCD; pixel 
scale 0.304"; field of view 10.4' x 10.4') and through a stan¬ 
dard Bessell R filter. The point spread function (PSF) was in¬ 
tentionally defocused to a radius of ~ 20 physical pixels to 
minimize flat-field residual errors and avoid detector satura¬ 
tion. The exposure time was set to 90 s, resulting in a net 
cadence of —115 s when considering the overheads. Science 
frames were bias- and flat-field-corrected by standard proce¬ 
dures. Photom etric measurements were made with the STAR- 
SAT pipeline (iNascimbeni et al.l 1201 lL 120131) . STARS KY de¬ 
livered the best differential light curve of KELT -6 (i.e. the one 
with the least scatter) using a set of four stable comparison 
stars (UCAC4 604-049448, UCAC4 604-049449, UCAC4 604- 
049450, and UCAC4 604-049454). Unfortunately, our obser¬ 
vations were plagued by technical problems that delayed the 
start time of the observations, and the telescope pointing was 
re-adjusted, causing an offset between the first and the second 
halves of the light curve, which was included as a parameter in 
the transit model. 


3. Stellar parameters 


The photospheric parameters were derived with different meth¬ 
ods from the co-added HARPS-N spectrum obtained from the 
out-of-transit observations (S/R~380 per pi xel at 5500 A). We 
used the LTE code MOOG (ISneden Il973l). along wi th atmo¬ 
spheric models dKuruczJ 19921 Castelli & Kuruczl2004l) and iron 
equivalent widths (EW). Two analyses were independently per¬ 
formed with the main differences being the list of iron lines 
and the technique used to mea sure EWs. In one case we 
used ARESv2 (ISousa et al.l 1201 5l) with th e automatic contin - 
i placement set-up and the line list f rom Sousa et akl ( 201T), 


adapted to the 2014 version of MOOG (Dumusaue etjd 


2014). 


Inthe other case, the EWs of the line list from Biazzo et al.l 
d2012l) were measured by hand using the IRAF task SPLOT, and 
the iron abundance was determined using the 2013 version of 
MOOG fsee iDamasso et akl 12015 ). In both cases, the analysis 
was performed differentially with respect to the Sun spectrum. 
We also did a third independent analysis based on the method de¬ 
scribed in lGandolfi et al.1 (|2015l) . We fitted th e HARPS-N spec¬ 
trum t o a grid of theoretical models from ICastelli & Kuruczl 
120041) . using spectral features that are sensitive to different pho¬ 
tospheric parameters. 

All the analyses gave consistent results, so we calculated 
their weighted averages and adopted the average of the in¬ 
dividual uncertainties as errors. Table Q] summarizes our re¬ 
sults. We notice that our best-fit value for T e $ is 170 K higher 
than value adopted by Col4, which was obtained from spec- 
tral synthesis mode lling with Spectroscopy Made Easy (SME; 

I Valenti & Piskunovll 19961) using HIRES spectra. This difference 
could arise from the different analysis techniques rather than 
from the properties of the spectra, by noting that a well known 


Table 1. Stellar parameters for the star KELT -6 derived from the analy¬ 
sis of the HARPS-N spectra and from stellar evolutionary tracks. Esti¬ 
mates from the KELT -6 b discovery paper are also listed. 


Parameter This work Co 14 Note 


T eff [K] 
log g [cgs] 

[Fe/H] [dex] 

Microturb. £ [km s -1 ] 
Vsin/* [km s _1 ] 

[i linear law limb-darkening coeff. 
Mass [M 0 ] 

Radius [R 0 ] 

Age [Gyr] 

Luminosity [L 0 ] 

Density [g cm -3 ] 


6272±61 

6102±43 


4.12±0.07 

4 074+ 0 045 
' -0.070 


-0.27±0.06 

0 981+0 039 
-U.Z8I_ 0038 


1.49±0.1 

0.85 (fixed) 


4.53±0.26 

5.0±0.5 

(i) 

0.48±0.14 


(i) 

1.126±0.058 

1 DR S+0.043 
l.U85_ 0040 

(2) 

1 S9Q+0143 
-0.137 

1 ^80+0160 

1 . joU_o 094 

(2) 

4 QO+0.66 

-0.46 

6.1±0.2 

(2) 

3.24±0.62 

oii +0.68 

-0.39 

(2) 

0 44+ 0 - 15 
u -^-o.io 

D T 87+0 068 

U. jo '-0.088 



Notes. (1) Derived from HARPS-N spectra (RM effect). (2) Matching T c gr, [Fe/H], and 
log g to the Yonsei-Yale evolutionary tracks. 


bias exists in some versions of SME (( Torres et al.l 20121). which 
was removed in more recent versions ( Brewer et''al]|2015i ). 

The stellar mass, radius, and age were determined by com¬ 
paring our measure d 7Vff, log q, and [Fe/Hj with the Yonsei-Yale 
evolutionary tracks ( Dernaruue et al.l2004l) through the ^ 2 3 statis¬ 
tics (ISanterne et al.l 1201 ll) . Results are listed in Table Q] The 
adopted errors include an extra 5% in mass and 3% in radius 
added in quadrature to the formal errors to take systematic uncer¬ 
tainties into account in the stellar models ([South worthll201 ll ). We 
also employed the stellar density p* derived by Co 14 as a p roxy 
of the stellar luminosity instead of log <7 (e.g. ISozzetti et al.l 
[2007 ). The results are fully consistent with the previous find¬ 
ings and not more precise. We thus adopted the parameters ob¬ 
tained using log g, finding that KELT -6 appears to be slightly less 
evolved than stated by Col4. 


4. Improved transit ephemeris for KELT-6b 

The light curve observed with the IAC-80 telescope is shown in 
the upper panel of Fig. Q] It has an average photometric scatter 
of 0.9 mmag on a 115 s timescale. 

Since we could not measure the out-of-transit flux of KELT- 
6 to properly normalize the off-transit level of our light curve, a 
full, detailed modelling to extract all the orbital parameter^ of 
KELT -6 b was hampered by the unsolvable degeneracy between 
the normalization level and the other transit parameters. There¬ 
fore we chose to fix R*/a, R p /a, and i to the best-fit values pub¬ 
lished by Co 14, limiting the number of free parameters to foufl 
The fit was performed with the code JKTEBOP v34 dSouthworthl 
120081 ) . and the associated errors were derived through 10000 
classical Monte Carlo iterations. 

From the fit we could only determine the time of central 
transit 7' c of KELT -6 b. We then derived an improved ephemeris 
(T c and orbital period P) from the observed-minus-calculated 
diagram (lower panel of Fig. [Tji. by fitting both our measure¬ 
ment of r c and the high-quality (“primary”) points from Col4 
through a weighted least-squares procedure, as well as an ad¬ 
ditional ephemeris obtained by the KELT team from an unpub- 

2 The fractional radii R^/a and R p /a, with R* being the star radius, R p 
the planet radius, and a the planet semi-major axis; the orbital inclina¬ 
tion 7 ; the orbital period P\ the time of central transit T 0 . 

3 i.e. T 0 \ the linear term u 1 of a quadratic limb darkening law, with the 
quadratic term interpolated from the lClaret & Bloemen! d201 ll) tables; 
the off-transit normalization level; an existing zero-point offset between 
the first and the second halves of the light curve. 
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Fig. 1. Upper panel: KELT-6 b transit light curve (2014 Apr 11. IAC-80 
telescope), and residuals from the best fit model (which is plotted with 
a solid blue line). The root mean square of the residuals is 0.9 mmag. 
Lower panel: O - C diagram including the highest-quality measure¬ 
ments from COM plus one unpublished measurement from the KELT 
team (blue points on the left side) and from our IAC-80 observation (red 
point). The O-C diagram is folded on our newly determined ephemeris 
(solid red line), while the Col4 ephemeris is plotted with a dotted blue 
line. 


lished transit that occurred on 2014 February 12. The results are 
shown in Table[4] with the reference epoch set on the most recent 
transit. 


5. Analysis of the radial velocities 

5.1. KELT-6c comes out 

The HARPS-N and TRES RVs collected for this work are listed 
in Tables [2] and [3] respectively. In our sample, we also in¬ 
cluded the RVs used in Co 14 and collected by TRES and the 
High Resolution Echelle Spectrometer (HIRES) at the Keck 
telescope, while we excluded the first 26 HARPS-N measure¬ 
ments taken during the night of the KELT-6 b transit. The RV 
time series shows a long-term modulation with a clear turnover 
and a semi-amplitude higher than that of the short-period sig¬ 
nal due to KELT-6 b (Fig. E}. KELT-6 appears to be a quiet star 
(<log(R^ K )>=-4.992, cr lo g (R ^ K )=0.021), without evidence in the 
log(R' ) data for a modulation ascribable to an activity cycle 
(Fig. [3J upper plot). Besides the log(RJ |K ) index, we calculated 
the CCF bisector span (BIS) from the HARPS-N spectra, and 
we found that no significant correlations exist over the times- 
pan of the HARPS-N measurements between these datasets and 
the RVs with the orbital solution of KELT-6 b removed (Fig. [3 
middle and lower plots). This evidence supports the hypothe¬ 
sis that the observed RV long-term variations are not due to 
the stellar activity but to an outer companion, which we name 
KELT-6 c and for which the data cover almost one complete or¬ 
bit. This scenario is further strengthened by looking, for exam¬ 
ple, at the empirical calibrations of lSantos et al.l rt2000l) (see Eq. 
2 and Fig. 6 therein) that, for a star like KELT-6, predict a RV 
dispersion due to activity-related phenomena at the level of ~10 
m s -1 , while our measured semi-amplitude is -3-5 times higher 
(K ~66 m s 1 ) once the orbital solution of KELT-6 b is removed. 



6000 6200 6400 6600 6800 7000 7200 

i . ; 

1 • ** > 

l Si .p* 


B 6000 6200 


6400 6600 6800 7000 

Time [BJD TDB-2,450,000] 


Fig. 2. Upper panel: Radial velocity time series of KELT-6. The l-<x 
errorbars are over-plotted to each data point. When they are not visible, 
this means they are smaller than the symbol size. Over-plotted is our 
best two-planet Keplerian model (orange line), calculated from the best- 
fit orbital parameters of Table[4] Lower panel: Residuals of the best-fit 
two-planet model. 


Orbital parameters and uncertainties for the KELT planets were 
determined with a Bayesian differential evolution Markov chain 
Monte Carlo analysis (e.g. iDesidera et ai]l2014l) . We adopted a 
two-planet Keplerian model with sixteen free parameters, as¬ 
suming for KELT-6 b Gaussian priors on Tq and P based on the 
new ephemeris. Fitted and derived parameters are listed in Ta¬ 
ble [4] The best-fit model is shown in Fig. [2] and the residuals 
show a dispersion of ~17 m s -1 , which reduces to 9.5 m s'" 1 
when considering only the more precise HIRES and HARPS- 
N measurements. This fact agrees nicely with the level of the 
activity-related jitter predicted bv ISantos et al.l d2000i) . 

KELT-6c appears to be massive (M p sin/=3.71+0.21 Mj up ) 
and moves on a moderately eccentric orbit (significant at the 
~5cr level) with period P ~3.5 years. Residuals of the two-planet 
model do not reveal evidence of any trend of having a third com¬ 
panion (Fig. [2j- Our results suggest that KELT-6 b has a likely 
circular orbit, with the eccentri city e = 0.029 4 5 ! 1 compatible 
with zero within ~2.45cr dLucv & Sweenevlll971 ). A reanalysis 
of the data done by the KELT team corrected the eccentricity to 
the value e-0. 058+0.034, which is in accordance with our mea¬ 
surement. 


5.2. The Rossiter-McLaughlin effect for KELTSb 

The RV time series covering the transit of KELT-6 b was anal¬ 
ysed using the numeri cal model and the least- squares fitting 
algorit hm described in ICovino et al.1 (|2013l) and lEsno sito et al. 
( 2014 ). Four parameters were set freqj, while all other relevant 
parameters were kept fixed to the values derived from the spec¬ 
troscopic and photometric analyses. The best-fit values are re¬ 
ported in Tables|T]and[4] together with their uncertainties derived 


4 The central transit epoch T 0 , the orbital period P, the RV semi- 
amplitude K, ^fe cos u>, and \fe sin at of both planets KELT-6 b and 
KELT-6 c (e and to being the eccentricity and the argument of perias- 
tron); a jitter term and a RV offset for each of the different datasets. 

5 y: barycentric RV at mid-transit epoch; A: projected angle between 
the planetary orbital axis and the stellar spin axis; Esin/,: projected 
stellar rotational velocity; p: limb-darkening coefficient of a linear law. 
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Table 4. Orbital and physical parameters for KELT-6 b and KELT-6 c. 


Parameter 

KELT-6 b 

KELT-6 c 

Primary transit analysis 

T c [BJDtdb -2,450,000] 

7124.50954 ± 5.7 ■ 10~ 4 


P [days] 

7.845582 ± 7 ■ 10” 6 


(Rp/Rsf 

0 077613+ 0 0010 


b°” 

0 20 +014 
^ -0.13 


i a - c [deg] 

RV analysis 

88.81/J™ 


P d [days] 

7.8455821 ±7-10- 6 

1276 + |* 

K a [m s~>] 

41.8 ± 1.1 

65 7 +2 - 6 

e e 

0-029/°;™ 

0 2l+° o 39 

u ‘ z ' 1 — 0.036 

( y/e sin co) d 

-0.083/S 

-0.454 ± 0.042 

(cos cj) d 

0 126 +0035 

U ' A -0.046 

-0.011 ±0.071 

D e [deg] 

308/™ 2 

268.7 ± 8.8 

T/[BJDtdb -2,450,000] 

7124.50954 ± 5.8 ■ 10 -4 

7432/™ 

yTREs'' [ms _I ] 

46.3 ± 3.8 


rmREsGnis- 1 ] 

22.1 ±8.5 


THARPS-l/ [ms -1 ] 

1128.9 ±3.0 


(Jitter-pRps) d [mr 1 ] 

14.2 ±4.5 


C/itterniREs)1 m s 1 1 

*■133 


U''tt«taARPS— nH [“S' 1 ] 

2.2 ± 1.2 


(MpSiniW [M ju p] 

0.442 ±0.019 

3.71 ±0.21 

R P S [Rjup] 

1.18 ±0.11 


V [AU] 

0.080±0.001 

2.39 ±0.11 

7rm' [ms" 1 ] 

1190.1 ± 1.4 


A 1 [deg] 

-36±11 



Notes. (o) Co 14 ^ Impact parameter. ^ Inclination of the orbital plane . (d) Fitted. De¬ 
rived. ^ For KELT-6 b this corresponds to the real mass of the planet, derived by assuming 
the best estimate of the inclination angle of the orbital plane provided by Co 14. ^ De¬ 
rived from the transit parameter R p /R s in Co 14 and our value of R s through a Monte-Carlo 
analysis. ^ Derived from the third Kepler’s law through a Monte-Carlo analysis, using 
our updated values for the stellar mass and planetary orbital periods. ^ Fitted (Rossiter- 
McLaughlin effect). 



-0.015 -0.01 -0.005 0 0.005 0.01 0.015 0.02 0.025 

phase 


Fig. 4. HARPS-N RV time series spanning the transit of KELT-6 b oc¬ 
curred on 2015 April 11 (open blue circles) and our best-fit model of 
the Rossiter-McLaughlin effect (red solid line). 


by means of a bootstrapping method. The best-fit RM model is 
shown in Fig. [4] Results show that KELT-6 b moves in a pro¬ 
grade orbit and is slightly misaligned with respect to the stellar 
spin axis with a projected spin-orbit angle /t=-36° + 11°. The 
lack of enough RV data in the pre-transit phase means that a 
value of A closer to zero cannot be completely ruled out, but our 
solution comes with a lower x 1 than for the case of A fixed to 0 
(1.39 vs. 1.89). 


6. Discussion 

KELT-6 joins a small number of host stars with a transiting 
planet and a measured RM effect for which there is evidence of 
outer companions. To our knowledge they are HAT-P-13, HAT- 
P-17, and WASP-8 (iKnutson et al.l[2014l) . the last two hosting a 
transiting planet with the orbital period close to that of KELT-6b 
(respectively, P= 10.3 and 8.2 days) and with WASP-8 being mis¬ 


aligned (A-- 123 .3^440. Our results allow the description of the 
KELT-6 system architecture in some detail, and considerations 
can be made about its dynamical evolution and stability. 

We tested the coplanar high-eccentricity migration hypoth¬ 
esis (CHEM, IPetrovichfeO 15a ) as a possible mechanism to ex¬ 
plain the current orbital elements of KELT-6. CHEM explores 
the possibility that hot Jupiters are formed through secular grav¬ 
itational interactions with an outer planet on an initial eccentric 
orbit, and then circularized by tides with the host star at perias- 
tron. The mechanism describes the planetary migration as occur¬ 
ring on nearly the same plane where the planets formed, assum¬ 
ing a low mutual inclination of the orbital planes (< 20°), and it 
predicts that hot Jupiters should have distant and more massive 
companions with moderately high eccentricities (e ~ 0.2-0.5). 
CHEM reproduces the observed architecture of KELT-6, pro¬ 
vided that the inner planet, initially with e>0.5, started migra¬ 
tion inside ~ 1 au, and there was an initial moderate mutual in¬ 
clination of ~10°-20° to account for the measured angle A. This 
suggests that the current mutual inclination of the orbital planes 
could be lying in this range. We note that the system should be 
old enough for the orbit of KELT-6b to have been circularized. 
In fact, the circularization time with the current system parame¬ 
ters is about 0.45 Gyr, assumin g a m odified tidal quality factor 
Q' p = 10 5 for the planet dLainev et ali 2 0091) . where most of the 
tidal dissipation is expected to have occurred in the case of an 
initially eccentric orbit (see iJackson et al. l 12008 ). We speculate 
that the initial conditions required by CHEM to explain this sys¬ 
tem could be the result of a preceding phase of planet-planet 
scattering. A system with more than two planets in nearly cir¬ 
cular and coplanar orbits at > 1 au might have become unstable, 
losing planets by ejections and leading to eccentricity excitation 
and moderate change (b y a factor up to ~2) in the semi-major 
axis of the inner planet (lJuric & Tremainel2008l:lChatteriee et al.l 
120081 ) . 

We investigated the dynamical stability of the KELT-6 sys¬ 
tem by using the empirical relation (17) in IPetrovichl ( 2015 b). 
who studied the final results of the evolution on long timescales 
of two-planet systems with arbitrary eccentricities and mutual 
inclinations against either ejections or collisions with the host 
star. We found that the relation is satisfied well, implying a sta¬ 
bility preserved on a timescale of a few Myr (i.e. the range of 
validity of Petrovich’s formula), as expected for a system that is 
already ~5 Gyr old. 
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Fig. 3. [Upperpanel). Time series of the log(RJj K ) chromospheric index 
as measured from HARPS-N spectra. [Middle and lower panels). Cor¬ 
relation analysis between the RV residuals, obtained by removing the 
Keplerian signal due to KELT-6 b from the original HARPS-N dataset, 
and two indicators of stellar activity derived from the HARPS-N spec¬ 
tra: the log(RJj K ) chromospheric index and the CCF bisector span. Over¬ 
plotted are the values of Spearman's rank correlation coefficients, which 
show the absence of significant correlations and support the hypothesis 
that the origin of the observed RV long-term variations is Keplerian. 


Table 2. Sample list of all the radial velocities (RV) measured for KELT- 
6 with the HARPS-N spectrograph in the framework of the GAPS pro¬ 
gramme. Also indicated are the bisector velocity span (BIS) and the 
activity index log(RJj K ). The BIS are those derived by the data reduc¬ 
tion pipeline of HARPS-N with uncertainties assumed to be twice those 
on the radial velocities. 


Time 

RV 

RV err. 

BIS 

1o §(*hk) 

l°g(/?' HK ) err. 

(BJDj/t-c - 2,450,000) 

(m s -1 ) 

(m s" 1 ) 

(m s -1 ) 



6698.624617 

1127.6 

3.7 

85.6 

-4.953 

0.017 

6699.598118 

1124.1 

2.7 

92.6 

-4.995 

0.011 

6701.572178 

1069.7 

3.2 

76.1 

-4.993 

0.016 

6702.578027 

1044.5 

3.6 

83.2 

-4.967 

0.015 

6728.651408 

1105.1 

3.8 

91.5 

-5.020 

0.021 

6786.609623 

1143.6 

4.6 

103.4 

-4.952 

0.027 

6787.608024 

1112.3 

3.9 

101.6 

-4.983 

0.023 

6817.482897 

1176.6 

2.7 

83.9 

-5.022 

0.013 

6819.529437 

1115.1 

3.6 

78.6 

-4.987 

0.019 

6820.503665 

1101.6 

5.3 

78.4 

-5.035 

0.040 

6861.396206 

1150.9 

2.5 

87.9 

-5.001 

0.011 

6874.367191 

1144.3 

3.7 

94.1 

-4.996 

0.020 

6875.369909 

1119.5 

9.1 

67.0 

-4.975 

0.076 

6987.775868 

1204.0 

3.6 

86.8 

-4.968 

0.019 

6998.737459 

1198.9 

5.5 

81.9 

-4.950 

0.032 

6999.772079 

1166.2 

5.7 

107.6 

-5.021 

0.044 

7005.756862 

1234.1 

9.3 

71.0 

-5.034 

0.086 

7011.773901 

1226.4 

2.7 

90.9 

-4.997 

0.012 

7028.700777 

1240.1 

3.5 

101.0 

-5.015 

0.019 

7029.691968 

1216.8 

3.5 

87.4 

-5.001 

0.018 

7050.675470 

1224.3 

4.8 

84.4 

-4.985 

0.030 

7068.780530 

1218.8 

4.7 

109.7 

-4.971 

0.028 

7109.622389 

1170.3 

5.0 

87.8 

-5.000 

0.031 

7112.604654 

1155.1 

8.0 

78.7 

-4.996 

0.069 

7113.520997 

1209.9 

3.7 

86.7 

-4.977 

0.021 

7114.443438 

1212.0 

5.0 

74.6 

-4.960 

0.032 

7115.451479 

1219.3 

2.6 

90.4 

-5.003 

0.012 

7124.384992 

1203.2 

4.1 

89.8 

-4.958 

0.024 

7124.395594 

1203.2 

3.7 

95.3 

-4.983 

0.022 

7124.406149 

1206.1 

3.3 

101.8 

-4.989 

0.018 

7124.416785 

1217.2 

3.2 

109.7 

-5.013 

0.019 

7124.427595 

1213.6 

3.4 

97.6 

-4.994 

0.020 

7124.438266 

1210.4 

3.3 

85.0 

-4.992 

0.018 

7124.448937 

1215.6 

3.4 

78.9 

-4.968 

0.018 

7124.459481 

1212.7 

3.4 

99.1 

-5.013 

0.020 

7124.470175 

1214.5 

3.6 

86.5 

-5.006 

0.022 

7124.481309 

1198.8 

3.0 

82.7 

-5.013 

0.016 

7124.491679 

1202.4 

2.9 

94.5 

-4.994 

0.014 

7124.502581 

1198.6 

3.1 

80.4 

-5.001 

0.016 

7124.513125 

1196.5 

3.2 

85.3 

-4.987 

0.016 

7124.524155 

1186.0 

3.1 

83.6 

-5.006 

0.016 

7124.534618 

1187.8 

3.3 

84.4 

-4.988 

0.018 

7124.543668 

1173.8 

7.7 

97.4 

-4.956 

0.056 

7124.556156 

1179.0 

3.2 

91.7 

-5.013 

0.018 

7124.566665 

1180.4 

3.1 

98.1 

-5.005 

0.017 

7124.577441 

1171.2 

3.1 

92.9 

-5.015 

0.017 

7124.587949 

1176.1 

3.2 

89.6 

-4.997 

0.017 

7124.599130 

1174.4 

3.3 

72.1 

-5.001 

0.019 

7124.609685 

1170.2 

3.2 

91.4 

-5.011 

0.018 

7124.620229 

1182.9 

3.1 

103.7 

-4.998 

0.017 

7124.631085 

1181.3 

3.4 

83.3 

-4.996 

0.018 

7124.641640 

1183.5 

3.4 

69.9 

-5.002 

0.018 

7124.652496 

1189.0 

3.4 

91.8 

-4.976 

0.017 

7124.663202 

1176.4 

3.5 

81.3 

-4.982 

0.018 

7124.673908 

1188.7 

3.3 

83.9 

-4.985 

0.016 

7124.684347 

1186.2 

3.2 

89.2 

-4.987 

0.017 

7124.695354 

1182.3 

2.9 

101.9 

-4.995 

0.015 

7124.705932 

1179.9 

2.9 

101.8 

-4.980 

0.015 

7142.537170 

1144.0 

3.6 

89.5 

-4.990 

0.021 

7143.470202 

1158.9 

3.4 

86.3 

-5.018 

0.020 

7144.467594 

1196.6 

3.8 

84.4 

-5.003 

0.023 

7150.499394 

1155.1 

4.6 

90.2 

-5.004 

0.026 

7159.414630 

1164.7 

3.2 

99.1 

-5.024 

0.019 

7170.469914 

1214.6 

2.7 

83.9 

-5.010 

0.012 

7173.479806 

1143.8 

5.2 

112.2 

-4.961 

0.033 

7174.499116 

1140.4 

5.3 

91.6 

-4.937 

0.032 

7175.462458 

1179.1 

2.7 

99.2 

-4.996 

0.013 

7176.476235 

1212.3 

3.8 

82.5 

-4.994 

0.021 

7205.414495 

1136.2 

2.5 

90.2 

-5.000 

0.011 

7206.421521 

1157.3 

3.1 

80.1 

-4.968 

0.014 

7207.432852 

1189.7 

3.8 

87.3 

-4.952 

0.021 
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Table 3. List of the new relative radial velocities of KELT-6 measured 
with the TRES spectrograph and used for the first time in this work. 


Time 

(BJD[/7 -c - 2,450,000) 
6640.021901 
6643.998845 
6698.953413 
6702.883586 
6729.914621 
6733.902128 
6801.709729 
6803.711379 
6804.800386 
6812.690753 
6816.685639 
6824.685293 
6993.029954 
7055.946535 
7090.798944 
7094.908058 
7110.916771 
7121.849369 
7141.862950 
7145.767451 
7149.696112 
7169.684005 


RV 

RV error 

(m s- 1 ) 

(m s’ 1 ) 

-30.50 

36.12 

35.66 

23.47 

86.48 

24.29 

-15.28 

28.33 

64.77 

24.93 

4.10 

19.52 

59.28 

22.98 

-15.48 

12.28 

-14.19 

22.87 

54.72 

29.16 

73.70 

18.88 

94.65 

14.15 

107.42 

22.14 

94.92 

18.22 

172.10 

19.99 

101.59 

19.32 

73.25 

20.13 

133.02 

17.07 

35.84 

17.49 

156.82 

17.82 

42.03 

18.32 

133.14 

18.14 



